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The papain-like protease of severe acute respiratory syndrome coronavirus (PLpro) (EC 
3.4.22.46) is essential for the viral life cycle and therefore represents an important antiviral 
target. We have identified 6MP and 6TG as reversible and slow-binding inhibitors of SARS- 
CoV PLpro, which is the first report about small molecule reversible inhibitors of PLpro. The 
inhibition mechanism was investigated by kinetic measurements and computer docking. 
Both compounds are competitive, selective, and reversible inhibitors of the PLpro with K is 
values ~10 to 20 |jlM. A structure-function relationship study has identified the thiocarbonyl 
moiety of 6MP or 6TG as the active pharmacophore essential for these inhibitions, which has 
not been reported before. The inhibition is selective because these compounds do not exert 
significant inhibitory effects against other cysteine proteases, including SARS-CoV 3CLpro 
and several cathepsins. Thus, our results present the first potential chemical leads against 
SARS-CoV PLpro, which might be used as lead compounds for further optimization to 
enhance their potency against SARS-GoV. Both 6MP and 6TG are still used extensively in 
clinics, especially for children with acute lymphoblastic or myeloblastic leukemia. In light of 
the possible inhibition against subset of cysteine proteases, our study has emphasized the 
importance to study in depth these drug actions in vivo. 

© 2008 Elsevier Inc. All rights reserved. 


1. Introduction 

Severe acute respiratory syndrome (SARS) is a life-threaten¬ 
ing atypical pneumonia caused by the SARS coronavirus 
(SARS-CoV) [1]. During the outbreak in 2003, this virus 
infected more than 8000 people, and the fatality rate in 
humans was as high as 15% (World Health Organization). Its 


easy transmissibility between humans makes the reemer¬ 
gence of SARS-CoV a distinct possibility, which entails an 
urgent need to prepare antiviral drugs against SARS-CoV. To 
date, there is no effective anti-SARS drug available with 
proven efficacy in vivo [2], although vaccines and mono¬ 
clonal antibodies against SARS-CoV are being developed 
[3,4]. 
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SARS-CoV is a plus-strand RNA virus encoding four 
structural proteins, 16 nonstructural proteins, and 8 accessory 
proteins [5]. The 16 nonstructural proteins are the cleavage 
products of the 2 large polypeptides ppla and pplab, generated 
by the virally encoded proteases 3-chymotrypsin-like protease 
(3CLpro) and papain-like protease (PLpro) (EC 3.4.22.46). SARS- 
CoV also encodes a helicase with multiple enzymatic activities 
in vitro, including RNA/DNA helicase, nuclease, and RNA 5'- 
triphosphatase activities, and an RNA-dependent RNA poly¬ 
merase [6]. Extensive and in-depth investigation of SARS-CoV 
has led to the identification of its cellular receptor, carbox- 
ypeptidase angiotensin-converting enzyme 2 [7]; the determi¬ 
nation of the structures of many viral structural and 
nonstructural proteins [5]; and an understanding of the 
functions of many of these viral proteins in vivo [6,8]. 

An intensive hunt for effective anti-SARS drugs has been 
undertaken, including the screening of existing antiviral or 
other marketed drugs for possible anti-SARS effects [9]. 
Inhibitors have been both designed and discovered by targeting 
different components of the virus, including the interaction 
between the viral S protein and the host receptor, which 
mediates fusion and viral entry [10], cathepsin L (to decrease 
membrane fusion) [11], 3CLpro [9], RNA-dependent RNA 
polymerase [12,13], and helicase [14,15]. So far, these com¬ 
pounds have been predominantly tested for their in vitro 
efficacy and no human clinical trials have ever been initiated. 
Recently, an animal model mimic the mortality and human 
disease syndrome as induced by a mouse-adapted SARS-CoV 
has been set up [16], which will provide the opportunity to test 
the effectiveness of the inhibitors for in vivo efficacy. Proteases 
remain one of the most prominent and effective drug targets in 
antiviral therapies, including those directed against human 
immunodeficiency virus, hepatitis C virus, and other viral 
infections [17,18]. Because proteolytic processing is essential for 
the generation of a functional replication complex, 3CLpro and 
PLpro are potentially effective targets for anti-SARS drugs. Most 
efforts reported until now have focused on the development of 
3CLpro inhibitors [9]. However, except our previous observation 
that the zinc ion inhibits PLpro of SARS-CoV [19], no inhibitor of 
this protease has yet been reported. 

PLpro of SARS-CoV is a papain-like cysteine protease with 
deubiquitinating activity. It is located on polypeptide la and 
self-cleaves at three sites on the polypeptide. Its catalytic triad 
(Cysl651-Hisl812-Aspl826) has a broad range of pH optima, 
characteristic of papain-like cysteine proteases [19,20]. More¬ 
over, PLpro has deubiquitinating activity in vitro, and is 
capable of hydrolyzing diubiquitin, polyubiquitin, and syn¬ 
thetic ubiquitin peptide substrates [20,21]. Recently, its crystal 
structure has been solved, and it shows great similarity to 
those of two other mammalian deubiquitinating enzymes 
(DUBs), USP14 and herpesvirus-associated uiquitin-specific 
protease (HAUSP) [22-24]. Despite low sequence identity 
(around 10%), these three proteins share similar tripartite 
architectures, comprising finger, palm, and thumb domains. 
There is also a ubiquitin-like domain at the amino terminus of 
PLpro, which might be responsible for its association with 26S 
proteasomes, as has been observed for USP14 and a yeast DUB 
called ubp6 [25]. PLpro also has de-ISGylation activity [21]. 
ISG15 is a ubiquitin-like protein involved in the posttransla- 
tional modification of proteins. The biological functions of the 


deubiquitination and de-ISGylation activities of PLpro and its 
in vivo substrates are unknown at present. 

Previously, we developed a fluorescence resonance energy 
transfer (FRET) assay to screen the inhibitors of PLpro in a 
high-throughput format [19]. Here, we report the discovery 
that 6-mercaptopurine (6MP) and 6-thioguanine (6TG) inhibit 
SARS-CoV PLpro. 6MP and 6TG are purine analogues effective 
in the treatment of patients with acute lymphoblastic or 
myeloblastic leukemia [26,27], Their structure-activity rela¬ 
tionship and inhibition mechanisms were investigated. To our 
knowledge, these are the first chemical compounds reported 
to be reversible inhibitors of PLpro. The implications of using 
6MP and 6TG in antiviral therapies are discussed. 


2. Materials and methods 

2.1. Materials 

Ni-NTA agarose was from Qiagen, Inc., Valencia, CA, USA. The 
compound library used was the Genesis Plus Collection 
purchased from MicroSource Discovery Systems, Inc., Gay- 
lordsville, CT, USA. This library contains many compounds 
already approved by the United States Food and Drug Admin¬ 
istration (FDA). 6-Thioguanine (6TG), hypoxanthine, 6-methyl- 
mercaptopurine, and N-ethylmaleimide (NEM) were from 
Sigma-Aldrich Corporation, Inc., Louis, ST, USA. 6-Mercapto- 
purine (6MP) was from Fluka, Inc., Buchs, SG, Switzerland, and 
2-amino-6-methyl-mercaptopurine was from Acros Organics, 
Thermo Fisher Scientific Inc., Waltham, MA, USA. Papain-like 
cysteine protease cathepsins B, L, K and S were prepared as 
described previously [28-30] whereas papain was from Sigma- 
Aldrich Corporation, Inc., and further purified to homogeneity 
[31]. All calpains tested in this study were active titrated with E- 
64 from Sigma-Aldrich Corporation, Inc. [32]. 

2.2. Purification of SARS-CoV PLpro and 3CLpro 

The PLpro expression plasmid was constructed by amplifying 
the PLpro cDNA by PCR with the primers 5'-AAAGGATCCGCT- 
GAACTCTCTAAATG-3' and 5'-AAACTCGAGCGACACAGGCTT- 
GATG-3', using the previously constructed plasmid pBacPAK8- 
PLpro as a template [19]. The cDNA fragment was released 
with BamHI and Xhol digestion before it was ligated into the 
pET-22b vector with a 6x His tag at the C terminus (Novagen). 
PLpro was expressed in Escherichia coli BL21 (DE3). The cells 
were induced with 1 mM isopropyl-l-thio-p-D-galactose (IPTG) 
and cultured at 16 °C overnight, before they were collected by 
centrifugation at 6000 x g at 4 °C for 10 min. The cell pellets 
were sonicated in buffer containing 4 mM imidazole, 300 mM 
sodium chloride, and 20 mM Tris-HCl (pH 7.9), before the 
suspension was cleared by centrifugation at 10,000 x g for 
25 min at 4 °C. The supernatant was incubated with a 50% Ni- 
NTA slurry (Qiagen) at 4 °C for 60 min and loaded into a 
column. The column was washed with buffer containing 
30 mM imidazole, 300 mM sodium chloride, and 20 mM Tris- 
HCl (pH 7.9). PLpro protein was eluted with buffer containing 
400 mM imidazole, 300 mM sodium chloride, and 20 mM Tris- 
HCl (pH 7.9). The protein was stored in 50 mM phosphate 
buffer (pH 7.4), after buffer exchange with an Amicon Ultra-4 
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centrifugal filter device (Millipore). Typical yield of the protein 
was 2 mg per liter of cell culture. Purification and measure¬ 
ment of the enzymatic activity of 3CLpro of SARS-CoV were 
performed as described previously [33]. 


centrations of Dabcyl-FRLKGG and APIKGV-Edans peptides in 
a 1:1 ratio, and the enzymatic activity was assayed. The 
inhibition data were best fitted to a competitive inhibition 
pattern according to Eq. (2) 


2.3. Inhibitor-screening platform and deubiquitination 

assay 


Vmax [S] 

(l + [I]/K is )K m + [S] 



Inhibitor screening was carried out as described previously 
[19]. For the deubiquitination assay, PLpro (40 nM) was 
incubated with the chemical compounds for 10 min before 
the substrate, ubiquitin-AMC (1 |xM), was added in 50 mM 
phosphate buffer (pH 6.8). The enzymatic activities were 
determined by monitoring the enhanced fluorescence emis¬ 
sion upon substrate cleavage at excitation and emission 
wavelengths of 380 and 436 nm, respectively, in a PerkinElmer 
LS 50B luminescence spectrometer (USA). 

2.4. Steady-state kinetic analysis 

Because it is slightly more soluble in the assay buffer, we used 
Dabcyl-FRLKGGAPIKGV-Edans instead of Abz-FRLKGGA- 
PIKGV-Edans [19] as the substrate to measure the enzymatic 
activity of PLpro throughout the course of the study as described 
[19]. Specifically, the enhanced fluorescence emission upon 
substrate cleavage was monitored at excitation and emission 
wavelengths of 329 and 520 nm, respectively, in a PerkinElmer 
LS 50B luminescence spectrometer. Fluorescence intensity was 
converted to the amount of hydrolyzed substrate using a 
standard curve drawn from the fluorescence measurements of 
well-defined concentrations of Dabcyl-FRLKGG and APIKGV- 
Edans peptides in a 1:1 ratio. This will also correct for the inner 
filter effect of the substrate. For the kinetic analysis, the reaction 
mixture contained 0.5-25 |xM peptide substrate in 50 mM 
phosphate buffer (pH 6.2, 6.8, or 8.0) in a total volume of 
1 mL. After the addition of the enzyme to the reaction mixture, 
the increase in fluorescence at 520 nm (excited at 329 nm) was 
continuously monitored at 30 °C with a PerkinElmer LS 50B 
luminescence spectrometer. The increase in fluorescence was 
linear for at least 10 min, and thus the slope of the line 
represented the initial velocity (u). The steady-state kinetic 
parameters of the enzyme were determined by fitting the 
Michaelis-Menten equation (Eq. (1)) to the initial velocity data 


Vmax [S] 

K m + [S] 



in which V max is the maximum reaction rate, [S] denotes the 
substrate concentration, and K m is the Michaelis-Menten con¬ 
stant for the interaction between the peptide substrate and the 
enzyme. 

For the inhibition studies, 5-50 |xM peptide substrate and 0- 
100 [xM 6MP or 6TG were mixed with PLpro (0.1 |xM) at 30 °C, 
and the fluorescence intensity was monitored. Specifically, the 
enhanced fluorescence emission upon substrate cleavage was 
monitored at excitation and emission wavelengths of 329 and 
520 nm, respectively, in a PerkinElmer LS 50B luminescence 
spectrometer. Fluorescence intensity was converted to the 
amount of hydrolyzed substrate using a standard curve drawn 
from the fluorescence measurements of well-defined con- 


in which [I] is the inhibitor concentration and K is is the slope 
inhibition constant for the enzyme-inhibitor complex. 

2.5. Inactivation mechanism 

For the inactivation studies using 6MP, PLpro enzyme (0.1 |xM) 
in 50 mM phosphate buffer (pH 6.2) was incubated with 
different concentrations of 6MP (2-200 |xM) at 30 °C for 2 h. 
Aliquots were withdrawn at different time points and added to 
an assay mixture containing 50 |xM peptide substrate to 
determine the remaining enzyme activity. Substrate hydro¬ 
lysis in the presence of inhibitor was linear for the first 10 min, 
and a good fit of the pseudo-first-order rate equation (Eq. (3)) to 
experimental data was obtained 

-ln(E t /Eo) = feinactt (3) 

in which E 0 is the enzyme activity at time zero and E t is the 
enzyme activity at time t. The slope of this semi-log plot 
represents the observed inactivation rate constant (femact)- 
To study its inactivation by 6TG, PLpro enzyme (0.1 |xM in 
50 mM phosphate buffer (pH 6.2) was incubated with different 
concentrations of 6TG (5-100 |xM) at 30 °C) in the presence of 
the peptide substrate, and the enzymatic activity was traced 
for 10 min. All the progress curves recorded showed an expo¬ 
nential approach to a final linear rate and were analyzed by the 
least-squares fitting of the following integrated rate equation 
(Eq. (4); Ref. [34]) to the experimental data: 

[P] = u s t + U ' +Vs [1 - exp(-fe inact t)] + d (4) 

^inact 


in which Uj is the initial velocity, u s is the steady-state 
velocity, and d is the displacement on the y-axis. 

The replot of k inact versus 6MP or 6TG concentration was 
not linear and was therefore fitted to a saturation curve, 
according to Eq. (5) [34,35] 


"inact 


kmax [I] 
fQnact [1] 



in which K inact is the dissociation constant of the enzyme-6MP 
or enzyme-6TG complex and fe max is the maximum inactiva¬ 
tion rate constant. 


2.6. Computer modeling of PLpro in a complex with 6 MP or 

6TG 


The crystal structures of SARS-CoV PLpro (pdb code: 2FE8) [24] 
and 3CLpro (pdb code: 1UK2) were used as the templates. 
Docking was performed with DS Modeling 1.7 software 

O Q 

(Accelrys). Before docking, the size space of 110 A (891 grid 
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points) near the active site was chosen. The structures of 6MP 
and 6TG were created and separately specified as input 
ligands. During the docking experiment, the chemical flex¬ 
ibility of 6MP or 6TG was allowed, and the interaction energy 
was adopted as the force field and charging method. The grid 

o 

extension from the site was 3.0 A, and the nonbonded cutoff 
distance was 10.0 A. The softened potential energy and 
distance-dependent dielectric were used (the dielectric con¬ 
stant was 1.0). The docking results were calibrated by Monte 
Carlo trials, in which maximal torsion of 25 and possible poses 
of 10 were compared. 


3. Results 

3.1. Identification of inhibitors of SARS-CoV PLpro 

We expressed and purified PLpro from E. coli as described in 
Section 2. The purified protease ran at around 52 kDa on 
sodium dodecyl sulfate-polyacrylamide gel electrophoresis, 
with over 90% purity (data not shown). MASS analysis was also 
carried out to determine the homogeneity of the enzyme 
preparation. Single peak around the molecular mass of 52,346 
was observed, close to the predicted 52,645 Da (data not 
shown). When assayed with the peptide substrate under 
conditions similar to those reported previously [19], the PLpro 
from this preparation exhibited very similar kinetic properties 
to the PLpro isolated from insect cells with the specific activity 


of 130 [xMol min' 1 mg -1 using the substrate Dabcyl-FRLKGGA- 
PIKGV-Edans, as compared to that of 404 |xMol min' 1 mg' 1 
from insect cells using the substrate Abz-FRLKGGAPIKGV- 
Edans [19], suggesting that the expression system had no 
major influence on the properties. 

To identify potent inhibitors of SARS-CoV PLpro, we 
screened a library containing 960 compounds using a screen¬ 
ing platform that we had established previously [19]. Inter¬ 
estingly, we found that, apart from the zinc ion [19], 
thiocarbonyl-containing 6MP (compound 1) and 6TG (com¬ 
pound 2) were effective inhibitors of SARS-CoV PLpro, with 
IC 50 values of 21.6 and 5 |xM, respectively (Table 1). NEM 
(compound 6), a commonly used cysteine protease inhibitor 
that acts by covalently modifying the active-site Cys through 
Michael addition, was also found to be an effective inhibitor of 
PLpro with an IC 50 value of 4.4 |xM, however, it was not further 
investigated. Because PLpro is a deubiquitination enzyme with 
a structure highly homologous to those of other DUBs 
[20,21,24], we next tested whether 6MP and 6TG also inhibited 
the deubiquitinating activity of PLpro. Using ubiquitin-AMC as 
the substrate, we found that 6MP, 6TG, and Zn 2+ effectively 
inhibited the deubiquitinating activity of PLpro (Fig. 1). To 
determine whether this inhibition is reversible, we first 
incubated the enzyme with both 6MP and 6TG until no activity 
of the enzyme could be detected. Next, the enzyme was 
purified with gel filtration, followed by new activity measure¬ 
ment. The enzymatic activity was restored, suggesting that 
the enzyme-inhibitor complex dissociated on the column 
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Fig. 1 - Inhibition of PLpro deubiquitination activity by 
thiopurines and Zn 2+ . The percentages of PLpro activity 
remaining after treatment with 6MP, 6TG, or Zn 2+ were 
measured. The fluorogenic substrate ubiquitin-AMC 
(1 pM) was used as the substrate. 


and/or by dilution in the assay buffer, which is indicative of 
reversible inhibition. 

To identify the active pharmacophore in the inhibition, a 
structure-function relationship study was carried out. A sub¬ 
structure search based on the structures of 6MP and 6TG was 
carried out using the software ISIS (http://www.isis.rl.ac.uk/ 
dataanalysis/) against the MDL Available Chemical Directory 
(MDL ACD). Since the major difference between adenosine and 
6MP/6TG lies in the thiocarbonyl group, we searched for the 
different substituted group at this site in the database. 
Compounds 3, 4 and 5 were discovered and available to 
purchase commercially (Section 2). Replacement of the 
thiocarbonyl of 6MP or 6TG with either hydroxyl (compound 
3) or a methylthio group (compounds 4 and 5) resulted in 
compounds devoid of inhibitory activity (Table 1), suggesting 
that the thiocarbonyl group of 6MP or 6TG is the active 
pharmacophore for the inhibition. To gain further insight into 
the inhibition mechanism, a docking experiment was per¬ 
formed. PLpro exists as a monomer in solution [19]. In the 
crystallization form, there are three monomers existing in the 
asymmetric unit [24]. The active sites for all three units were 
the same (data not shown). So we chose subunit one for the 
modeling study. Both 6MP and 6TG fit well into the active-site 
cavity of PLpro. The sulfur atom of 6MP or 6TG was juxtaposed 
closely with the y-S of Cysl651 at a distance of 3.4 A, 
suggesting possible formation of a hydrogen bond (Fig. 2). 
This is consistent with the active chemical properties of the 
thiocarbonyl group. The best docking scores for 6MP and 6TG 
binding to PLpro were 23.9 and 24.4, respectively. Moreover, 
the binding energies of 6MP and 6TG to PLpro were 8.2 and 
14.7 kcal/mol, respectively. 

3.2. Inhibition mechanism 

To understand the kinetic mechanism of the interaction of 
6MP or 6TG with the enzyme, the enzymatic activity was 
measured with a series of substrate concentrations and 



Fig. 2 - Docking and inhibition mechanism of 6MP and 6TG 
with PLpro. 6MP (colored blue) and 6TG (green) could be 
targeted to the active site of PLpro. The distances between 
the sulfur atom of the active-site Cysl651 and the sulfur 
atoms of 6MP or 6TG were all 3.4 A. 


various inhibitor concentrations. The inhibition data were 
globally fitted to all possible kinetic models (competitive, 
noncompetitive, uncompetitive and mixed competitive). As 
shown in Fig. 3, a simple competitive inhibition pattern best 
described the data for the following two reasons: (i) the 
double reciprocal plot with all lines intercepted the y-axis for 
both 6MP (Fig. 3A) and 6TG (Fig. 3B) and (ii), the slope versus 
inhibitor concentration replot was linear (Fig. 3, both insets). 
Therefore, consistent with the docking result (Fig. 2), both 
6MP and 6TG may compete with the peptide substrate for the 
active site via a competitive mechanism. The best fit of 
Eq. (2) to experimental data yielded K is values of 19 and 
13 fxM for 6MP and 6TG, respectively (Table 2). As summar¬ 
ized in Table 2, in the pH range from 6.2 to 8.0, which was 
shown previously to be optimal for the enzymatic activity of 
PLpro [19], the K m values for substrate hydrolysis in the 
absence of inhibitors ranged from 17 to 25 |xM, whereas the 
k cat values ranged from 0.06 to 0.13 s _1 . The Michaelis 
constant (K m ) for the substrate and inhibition constant for 
6MP and 6TG (K is or K inact ) were similar. 

To understand the kinetic mechanism of this inhibition, a 
time-dependent change in enzyme activity as a function of the 
inhibitor concentration was determined by preincubating the 
PLpro enzyme with different concentrations of the inhibitors 
at 30 °C prior to measuring the enzyme activity. For 6MP, the 
semi-log plots of E t /E 0 versus time were linear in the initial 
stage (Fig. 4A), and the pseudo-first-order rate constant (kinact) 
at different inhibitor concentrations followed saturation 
kinetics (Fig. 4C). For 6TG, the emission curve was curved 
downwards, indicating a slow-binding mechanism (Fig. 4B). By 
fitting the data to Eq. (4) for 6TG (Fig. 4B) (Section 2), different 
kinact values at various concentration of 6TG were determined 
(Fig. 4D). When we further fitted the saturated curvature to 
Eq. (5), we determined that the K inact values for 6MP and 6TG 
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Fig. 3 - Competitive inhibition of PLpro by 6MP or 6TG. (A) 
6MP; (B) 6TG. The inset figures represent [6MP] or [6TG] vs. 
slope replots. The enzymatic activity of PLpro was 
measured in the presence of different peptide substrate 
concentrations (5-50 pM) and various inhibitor 
concentrations: 0 (closed circles), 2 (open circles), 10 (closed 
triangles), 25 (open triangles), 35 (closed squares), 50 (open 
squares), 100 (closed diamonds) pM. The solid lines are the 
best fit by the competitive inhibition equation (Eq. (2)). R sqr = 
0.920 and 0.929, respectively. The same experiments were 
repeated at least three times. The kinetic parameters such 
as K m , fe cat and K is from the best fit were shown in Table 2. 


were 48 and 32 |xM, respectively, and the fe max of inactivation 
were 0.12 and 0.92 ms -1 , respectively (Table 2). 6TG is a 
slightly more potent inhibitor than 6MP, as indicated by a 
slightly lower K inact and a eight-fold higher k max . In 
summary, both 6MP and 6TG are slow-binding inhibitors 
of SARS-CoV PLpro, although both 6MP and 6TG have binding 
affinities for PLpro similar to that of the peptide substrate 
(Table 2). Based on the above kinetic and modeling results, 
we propose the following kinetic mechanism (Fig. 5). 6MP or 
6TG might form a hydrogen bond with the enzyme’s active 
site Cysl651 through the thio moiety of the thiocarbonyl 
group. The hydrogen bonding with Cysl651 blocks the 
essential sulfhydryl group, preventing its acylation and 
inhibiting the enzyme’s activity. 

3.3. Activity of the inhibitors against other cysteine 
proteases 

Next, we investigated whether 6MP and 6TG inhibit other 
cysteine proteases. Initially, the other SARS-CoV protease 
3CLpro, which is also a cysteine protease, was tested for 
possible inhibition by 6MP and 6TG. A comparison with PLpro 
was therefore made at 200 |jlM inhibitor concentration and a 
2 h incubation time. Whereas both inhibitors at this concen¬ 
tration almost completely abolished the activity of PLpro (<5% 
residual activity), they were found to be much less effective 
inhibitors of 3CLpro (49% residual activity for 6MP and 25% for 
6TG, respectively). In the next step, we also used molecular 
modeling to dock the inhibitors to 3CLpro. The best docking 
scores for 6MP and 6TG to 3CLpro were 17.8 and 18.4, 
respectively, which are lower than those for 6MP and 6TG 
(23.9 and 24.4, respectively). Moreover, the binding energies of 
6MP and 6TG to 3CLpro were -14.7 and -16.2 kcal/mol, 
respectively, which are lower than those to PLpro (8.2 and 
14.7 kcal/mol, respectively). Finally, the sulfur atoms of 6MP 
and 6TG occur at distances of 5.4 and 8.5 A, respectively, from 
the 7 -S of Cysl45 of 3CLpro, too far away to maintain a 
functional hydrogen bond. Therefore, the lower binding 
energies, lower docking scores, and greater distances of the 
drugs from 3CLpro are consistent with the biochemical data 
that these two drugs are not effective inhibitors of 3CLpro. 
Moreover, 6MP and 6TG were not effective inhibitors of other 
cysteine proteases, including cathepsins B, K, L, and S, and 
papain, even at a concentration of 1 mM (data not shown). 
Therefore, it seems that 6MP and 6TG are reasonably selective 
inhibitors of PLpro and do not inhibit other cysteine proteases. 
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Fig. 4 - Time-dependent inactivation of PLpro by 6MP and 6TG. (A) Different concentrations of 6MP (2 pM, closed circles; 
25 pM, triangles; 35 pM, diamonds; 50 pM, squares; 100 pM, invert triangles; and 200 pM, open circles) was incubated with 
PLpro for 2 h. E t /E 0 indicates the residual enzymatic activity at time t. The solid lines show the best semi-log fits of the data 
at different 6MP concentrations by pseudo-first-order rate equation (Eq. (3)). (B) Different concentrations of 6TG (0 pM, 
closed circles; 5-100 pM, open circles) were incubated with PLpro for 10 min. The solid lines were the best-fit results 
according to the slow-binding equation (Eq. (4)). (C and D) The observed inactivation rate constants (feinact) from panels A and 
B were replotted against 6MP or 6TG concentration. The solid line is the fit of the saturation equation (Eq. (5)) (R sqr = 0.992 
and 0.975, respectively). The best-fit parameters (K inact and fe max ) are shown in Table 2. The same experiments were 
repeated at least three times. 
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Fig. 5 - Kinetic mechanism of inhibitors. The upper part 
denotes enzyme inhibition and inactivation by 6TG (or 
6MP). The lower part shows the putative hydrolysis of 
peptide substrate pathway. The original peptide substrate 
Dabcyl-FRLKGGAPIKGV-Edans (or ubiquitin-AMC) was 
cleaved at the Gly-Ala (or Gly-AMC) peptide bond. The N- 
terminal half (or AMC) was released as the first product 
while the C-terminal half (or ubiquitin) acylated the active 
site Cysl651, which was then deacylated and completed 
the catalytic cycle. 


4. Discussion 

6MP and 6TG are used clinically to treat leukemia. 6MP is a 
purine analogue, effective as an anticancer therapy because it 
is converted to 6TG nucleotides by hypoxanthine phosphor- 
ibosyl transferase in vivo. 6TG is then incorporated into DNA, 
preventing further replication. To our knowledge, 6MP and 
6TG are the first reversible chemical inhibitors documented 
to inhibit PLpro of SARS-CoV, whereas many inhibitors 
have been developed for the other protease of SARS-CoV, 
3CLpro [9], 

In this study, we have identified these two compounds as 
slow-binding inhibitors of SARS-CoV PLpro, likely compete for 
the same binding site as the substrate. At the current stage, we 
could not rule out the possibility that the inhibition is through 
binding to an allosteric site thereby changing the conforma¬ 
tion of the active site. Nevertheless, the inhibition is 
competitive, reversible and selective against PLpro. Although 
the affinity of 6MP and 6TG against PLpro is in the |xM range, 
they could be used as lead compounds for further optimization 
to enhance their potency against SARS-CoV. There are many 
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such examples reported in the past of successful antiviral 
treatments with old drugs [36]. In light of the availability of 
several animal models [4], especially the one infected by a 
mouse-adapted SARS-CoV, which reproduces many aspects of 
the human disease including morbidity, mortality and 
pulmonary pathology [16]. It would be interesting to test 
whether these two drugs are effective in blocking SARS-CoV 
replication in these animal models. 

Importantly, our structure-activity study (Table 1) has 
identified the thiocarbonyl group of 6MP and 6TG as the 
moiety for the inhibition of PLpro, which was not known 
before for these two compounds. In light of the acute 
toxicities associated with these drugs and the active 
thiocarbonyl against subsets of cysteine proteases, the in 
vivo action of 6MP and 6TG might be more complex than we 
realize at this stage. It is well documented that 6MP causes 
some adverse drug reactions in a significant number of 
patients, which have been found to be caused by a genetic 
polymorphism of thiopurine S-methyltransferase (TPMT). 
The polymorphism produces a TPMT with much lower 
activity [37,38]. These patients are less able to metabolize 
6MP to the inactive 6-methyl mercaptopurine, resulting in the 
accumulation of 6TG and subsequent life-threatening mye- 
losuppression [37-39]. About 1 in 300 individuals carries a 
mutated TPMT and does not express fully functional TPMT 
[40]. For these patients, less than 10% of the normal dosage of 
6MP can be tolerated [39] . It is likely that locally high dosage of 
6MP/6TG might inhibit the subset of cysteine proteases, 
thereby contributing to the severe cytotoxic effects observed 
by these two drugs. Furthermore, 6MP can also be metabo- 
lically processed by hepatic cytochrome P450 to purine-6- 
sulfenic acid, which binds to proteins by forming disulfide 
bonds with the thio moieties of the proteins both in vitro and 
in vivo [41,42]. Whether purine-6-sulfenic acid forms the 
disulfide bonds with cellular cysteine proteases remains to be 
studied. 

We speculate that DUBs are the in vivo targets of 6MP 
and 6TG, based on our study of the mechanism of inhibition 
of PLpro, the selectivity of 6MP and 6TG, and the structural 
similarity of PLpro to other DUBs. To date, there is no potent 
selective chemical compound available for the inhibition of 
cellular deubiquitination enzymes other than a ubiquitin 
aldehyde developed for laboratory studies of the function of 
DUBs by covalent attachment to the active-site Cys [43,44]. 
Hence, the biological consequences of inhibiting DUBs are 
not clear at the moment. Moreover, the functions of the 
deubiquitination and de-ISGylation activities of PLpro in 
vivo are also unknown, as are those of its in vivo cellular 
substrates. Adenovirus encodes a deubiquitinating enzyme 
[45], which is necessary for processing viral precursor 
proteins during virion maturation [46]. ISGylation and 
de-ISGylation are involved in the innate immune response 
to viral infection [47]. Therefore, whether 6MP and 6TG 
inactivate the cellular DUBs that contribute to adverse 
drug reactions or other biological consequences requires 
further investigation. Because both drugs are still used 
extensively in clinics, especially for children with acute 
lymphoblastic or myeloblastic leukemia [26,27], our study 
has emphasized the importance to study in depth these 
drug actions in vivo. 
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